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Abstract. There is no direct evidence for radiation domination prior to big-bang nucleosynthesis, 
and so it is useful to consider how constraints to thermally-produced axions change in non-standard 
thermal histories. In the low-temperature-reheating scenario, radiation domination begins as late as 
^ 1 Me V, and is preceded by significant entropy generation. Axion abundances are then suppressed, 
and cosmological limits to axions are significantly loosened. In a kination scenario, a more modest 
change to axion constraints occurs. Future possible constraints to axions and low-temperature 
reheating are discussed. 
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INTRODUCTION 

If the axion has mass > 10~^ eV, it will be produced thermally, with cosmolog- 
ical abundance Q.^h^ = (ma/130 eV) (10/^*f) , where is the effective number of 
relativistic degrees of freedom when axions freeze out [1-7]. If ma ^ 1 eV, axions free- 
stream to erase density perturbations while they are relativistic, and thus suppress the 
matter power spectrum on small scales, much like neutrinos [4, 6, 8-17]. Data from 
large-scale structure (LSS) surveys and cosmic microwave-background (CMB) observa- 
tions impose the constraint ma ^ 1 eV to light hadronic axions [18-20]. We restrict our 
attention to hadronic axions in this work. 

These constraints were determined in the standard radiation-dominated scenario. The 
transition to radiation domination after inflation might be gradual [21]. In a modified 
thermal history, relic abundances may change, due to modified freeze-out temperatures 
and suppression from entropy generation. 

The universe could have reheated to a temperature as low as 1 MeV [22-26]. This 
low-temperature reheating (LTR) scenario may be modeled simply through the decay of 
a massive particle ^ into radiation, with fixed rate . This decay softens the scaling 
of temperature T with cosmological scale factor a, increasing the Hubble parameter 
H{T) and leading to earlier freeze-out for certain relics. Entropy generation then highly 
suppresses these relic abundances. Kination models offer another alternative to the 
standard thermal history, without entropy production, and cause more modest changes 
in abundances [27]. 

Past work has determined the relaxation in constraints to neutrinos, weakly interacting 
massive particles, and non-thermally produced axions are relaxed in LTR [24, 28, 29]. 
Here, we present new constraints to thermally-produced axions in the LTR scenario. 
We point the reader to Ref. [30] for a discussion of the more modest changes to axion 



constraints in the kination scenario, and for additional details relevant to the following 
discussion. We conclude by discussing the impact of future LSS surveys and CMB 
measurements of the primordial helium abundance on the allowed parameter space for 
axions. 



LOW-TEMPERATURE REHEATING (LTR) 

In the LTR scenario, the density of ^ particles and radiation obey [24, 28, 31]: 

where p^j, and pR denote the energy densities in the scalar field and radiation, and a is 
the cosmological scale factor, whose evolution is given by the Friedmann equation. The 
reheating temperature Tj-h is defined by = i/4;r3g^rh/45 T^/Mp[ [6, 28, 31], where 
Mpi is the Planck mass and g^rh is the value of g* when T = T^h. 

At the beginning of reheating, (j) dominates the energy density. The temperature is 
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related to the radiation energy density by T p^' [6] . We integrate Eqs. (1) to obtain 
the dependence of T on a [30]. When the scalar begins to decay, the temperature rises 
quickly to a maximum and then falls as T oc a~^/^. This shallow scaling of temperature 
with scale factor results from the transfer of scalar-field energy into radiation. When 
Pr overtakes p^ near T ~ Trh, the epoch of radiation domination begins, with the usual 
T oc scaling [30]. 

During reheating, H{T) oc {T / T,^f / M^x [24, 28], the universe thus expands faster 
than during radiation domination, and the equilibrium condition T = n (ov) > H is 
harder to meet. Relics with freeze-out temperature Tp > Tmax will thus have suppressed 
abundances because they never come into chemical equilibrium. Relics with T^h ^ 7f < 
Tmax come into chemical equilibrium, but their abundances are reduced by entropy 
production. 



AXION PRODUCTION 

Standard hadronic axions with > 10~^ eV are produced by the channels + n~ ^ 
a + 7t^, 7t^ + ^ n^ + a, and 7t^ + n:^ ^a + 7t^[4, 6, 18, 32]. Numerically evaluating 
the expression from Ref. [4] for the axion-production rate T and solving Eq. (1) for 
H{T), we estimate the axion freeze-out temperature Tp using the condition r{TF) ~ 
H{Tp). As is lowered, axions freeze out earlier due to the higher value of H, as 
shown in Fig. 1 . As Tj-h increases, the T oc a^^/s epoch becomes less relevant, and Tp 
asymptotes to its standard value. Now, since F oc f-^ oc ml [4], higher-mass axions keep 
up with the Hubble expansion for longer and generally decouple at lower temperatures. 
The resulting axion abundance is [30] 
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FIGURE 1. The left panel shows the freeze-out temperature of the reactions 7r+ + tt" <-> tt*^ + a, 
7r+ + TT^ <-» 7r+ + a, and n~ n~ + a, as a function of r^h, for 4 different axion masses. The right 

panel shows the axion abundance ila normalized by its standard value il^ for 4 different axion masses. 



where ma,ev is the axion mass in units of eV. 

When Tj-h ^ Tp, the present mass density in axions is severely suppressed, because of 
entropy generation. Using the numerical solution for a{T), we obtain fla- In the right 
panel of Fig. 1, we show D.^ normalized by its standard value, D.^, as a function of Tj-h- 
For » Tp, the axion abundance asymptotes to Q.^. 

CONSTRAINTS TO AXIONS 

Most constraints to the axion mass come from its two-photon coupling gayy [5, 6, 33- 
42]. This coupling depends on the up-down quark mass ratio r, for which there are 
experimentally allowed such that gayy vanishes, and so constraints to axions from star 
clusters, helioscope, RF cavity, and telescope searches may all be lifted [7, 43]. In 
contrast, the hadronic couplings do not vanish for any experimentally allowed r values. 
Axion searches based on these couplings are underway, and have already imposed the 
ma ^ 1 keV range [44]. These couplings also determine the relic abundance of axions, 
and so constraints may be obtained from cosmology. 

Mass constraints to thermal axions from cosmology are considerably relaxed because 
of entropy generation. A conservative constraint is obtained by requiring that axions 
not exceed the matter density of Q.mh^ ~ 0.135 [16] and is shown by the dot-dashed 
hashed region in Fig. 2. If ^ 40 MeV, constraints are considerably relaxed. When 
^ 95 MeV, we obtain ma ^ 22 eV, equal to the standard result. 

Axions will free stream at early times, decreasing the matter power spectrum on length 
scales smaller than the comoving free- streaming scale, evaluated at matter-radiation 
equality: 



Afs ~ (196 Mpc/ma,ev) (TjTy) { 1 + In [0.45ma,eV (Ty/T^)] } . 



(3) 
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FIGURE 2. The left panel shows upper hmits to ma in the LTR scenario. The dot-dashed hatched region 
shows the region excluded by the constraint Q.ah?' < 0.135. The solid hatched region shows the axion 
parameter space excluded by WMAPl/SDSS data. At low Trh, limits to the axion mass are loosened. 
The right panel shows the estimated sensitivity possible with LSST measurements of the matter power 
spectrum [48, 49], or from Lyman- a forest measurements of clustering on smaller length scales [50]. The 
hatched region indicates the region excluded by WMAPl/SDSS measurements. 



This suppression is given by AP/P ~ — S^a/^m if <S [18, 19, 45] and imposes 
a constraint to Q.i^h^. Including entropy generation, the relationship between the effective 
axion temperature and the neutrino temperature Ty is 
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if 7f > Tih, ^ ^ (10.75/g,F)'/^ if Tp < T^h- (4) 



Using Sloan Digital Sky Survey (SDSS) measurements of the galaxy power spectrum 
[46] and Wilkinson Microwave Anisotropy Probe (WMAP) [47] l^'-year measurements 
of the CMB angular power spectrum, Refs. [18, 19, 45] derived limits of ma ^ 1 eV. We 
map these results into the {flah^-, ■^s) plane. 

We calculate Q.^{T^-^,me)h^ and Xi^{T^\i,ms) for axions in LTR, and thus obtain the 
upper limit to the axion mass as a function of 7^^, shown in Fig. 2. For this data set, the 
smallest length scale for which the galaxy correlation function can be reliably probed is 
Amin = 40 h-^ Mpc [18, 19, 30]. For < 35 MeV, Xfs < 40 h'^ Mpc, and this axion 
mass constraint is lifted. At high Tfh, the constraint from LSS/CMB data {O-ah^ % 0.006) 
supercedes the constraint ^Ji^ < Q.mh^- 

Future instruments, such as the Large Synoptic Survey Telescope (LSST), will mea- 
sure the matter power-spectrum with unprecedented precision (AP/P ~ 0.01) [48, 49]. 
This order of magnitude improvement over past work [5 1 , 52] leads to the improved sen- 
sitivity shown by the dotted line in Fig. 2. To estimate possible constraints to axions from 
LSST measurements of the power spectrum, we recalculated our limits using the approx- 
imate scaling AP/P ~ —SQ-n/Slta, assuming AP/P ~ 10-2 for A > Amin = 40 h'^ Mpc. 

We also estimate the possible improvement offered by including information on 
smaller scales (Amin ~ 1 Z?^^ Mpc), as probed by measurements of the Lyman-ct flux 
power spectrum [50], also shown in Fig. 2. This is indicated by the dashed line in Fig. 2. 
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FIGURE 3. Total effective neutrino number for axions with masses saturating the tightest 

bound on axion masses from Fig. 2. The thick black line indicates the anticipated sensitivity of CMBPol 
[55] to Ny^, through measurements of Yp. 



We can see that higher and lower T± values are probed because of information on 
smaller length scales. 

AXIONS, LTR, AND BBN 

Future limits to axions may follow from constraints to the total density in relativistic 
particles at T ~ 1 MeV. This is parameterized by the axionic contribution to the total 
effective neutrino number A^^^^ [4, 30]: 



Nf = 3 + -i4] ^(rF/r,h),^(y) 
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ify> 1, 



1 [^*,F-1]"^/' ifj«l. 



(5) 



For sufficiently high masses, the axionic contribution saturates to 5Ny^^ = 4/7 at high 

r^h [4]. In Fig. 3, we show A^f (T^-h), 
best cosmological bound for a given T^h. 



Trh [4]. In Fig. 3, we show (T^h), evaluated at whichever nia which saturates the 



A comparison between the abundance of He (Yp) and the predicted abundance from 
BBN places constraints Ny^ at T ~ 1 MeV [53]; thus constraints to ^He abundances are 
also constraints on ma and T^^. Here we apply the scaling relation [54]: 

ANf = y I (6.25Ay^ + 1 )' - 1 } . (6) 

Direct measurements of Yp, including a determination of Q.h from CMB observations, 
lead to the 68% confidence level upper limit of Ny^ < 3.85 [56-58]. From Fig. 3, we 



see that this bound cannot constrain or Tj-h. If future measurements reduce systematic 
errors, constraints to will be obtained for the lighter-mass axions. 

Constraints to ma and Tj-^ may follow from CMB measurements of Yp. ^Ke affects 
CMB anisotropics by changing the ionization history of the universe [59]. CMBPol (a 
proposed future CMB polarization experiment) is expected to approach AYp = 0.0039, 
leading to the sensitivity limit A^^^^ < 3.30 [55, 58-61]. As shown in Fig. 3, for T^h ^ 
15 MeV, such measurements of Yp may impose stringent limits on the axion mass. Also, 
if axions with mass in the eV range are directly detected, Yp might impose an upper limit 
to Trh. 

CONCLUSIONS 

LTR suppresses the abundance of thermally-produced axions, once T^^ ~ 50 MeV, as 
a result of dramatic entropy production. The cosmologically allowed window for 
is extended as a result. Future probes of the matter power spectrum or the primordial 
helium abundance may definitively explore some of this parameter space. 
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